Extended RPA with ground-state correlations in a solvable model 
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The ground states and excited states of the Lipkin model hamiltonian are calculated using a 
new theoretical approach which has been derived from an extended time-dependent Hartree-Fock 
theory known as the time-dependent density-matrix theory (TDDM). TDDM enables us to calculate 
correlated ground states, and its small amplitude limit (STDDM), which is a version of extended 
RPA theories based on a correlated ground state, can be used to calculate excited states. It is found 

, that this TDDM plus STDDM approach gives much better results for both the ground states and 

' the excited states than the Hartree-Fock ground state plus RPA approach. 
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^jPj. The random phase approximation (RPA) based on the Hartree-Fock (HF) gromid state and the quasi-particle 
^ RPA (QRPA) based on the Hartree-Fock-Bogoliubov (HFB) ground state when pairing correlations are important 
are standard microscopic and self-consistent theories which have extensively been used to study nuclear collective 
excitations 1]. RPA and QRPA are regarded as mean-field theories because they can be formulated as the small 
t , amplitude limits of the time-dependent HF theory (TDHF) and the time-dependent HFB theory, respectively. On 
the other hand the time-dependent density-matrix theory (TDDM) 2] which we have developed in the past decade 
, incorporates the effects of two-body correlations into TDHF, and its small amplitude limit (STDDM) [2j becomes 
J> a natural extension of RPA which includes two-body amplitudes and also the effects of ground-state correlations: 
The correlated ground state is obtained as a stationary solution of the TDDM equations. We have recently applied 
this TDDMplus STDDM approach to study the ground states and low-lying quadrupole states in unstable oxygen 
isotopes Although those results are quite encouraging, the validity of the TDDM plus STDDM approach should 

' be tested in further applications. In this letter, we apply TDDM and STDDM to the Lipkin model and calculate the 
I ground states and excited states. We will show that TDDM and STDDM give much better results than HF and RPA. 
' The correlated ground state |<i>o) in TDDM is constructed such that the two equations, 

Fi{af3)^{<^o\[a+a0,H]\^o)^O, (1) 
3 ■ F2{al3a'P') - {'fo\Ka+ap,a^, , H]\^o) = 0, (2) 
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are satisfied, where [ ] is the commutation relation and H is the total hamiltonian consisting of the kinetic energy term 
and a two-body interaction. Eqs.(^ and JSJ describe the conditions that the occupation matrix — ($o|a^aQ I'i'o} 
^ I and the two-body correlation matrix C^^^i^, = {^o\a'^,a'lj,ai3aa\^o) ~ •^{'''^aa'^$f3') ' where A is the antisymmetrization 
^ operator, are time-independent. We choose the single-particle wavefunction ipa as an eigenstate of the mean-field 
hamiltonian ho{po) which is a functional of the one-body density matrix po. The expressions for Eqs.0J and lO are 
given in Ref.'i]. We have shown in Ref.[^ that Eqs.QJ and ^ can be solved using the iterative gradient method: 
Starting from the HF ground state where n^^ = 5/3a(0) for an occupied (unoccupied) state and C^^^,,^, — 0, we 
iterate 

n"{N + _ f n°{N)\ _ / a c VV Fi(7V) 
C"{N + 1) ) ^ [C"{N) ) "^[bdj [f2{N) 

until convergence is achieved. The matrix elements a, b, c, and d, which are functional derivatives of Fi and F2 and 
depend on the iteration step N, are equivalent to those appearing in the hamiltonian matrix of STDDM. They are 
also given in Ref.Q. We have to introduce a small parameter a to control the convergence process. 

Excited states are obtained solving the STDDM equations for the one-body transition amplitude Xap{p) = 
i^oW'fjO.al^n) and the two-body transition amplitude Xa/sa'p'ip) — {^oWa''^'f3'^P'^o:\^ti) , where 1$^) is the wave- 
function for an excited state with excitation energy w^. The equations in STDDM can be written in matrix form 
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When the ground-state |$o) is approximated by the HF one and only the particle (p) - hole (h) and 2p - 2h amplitudes 
(and their complex conjugates) are taken in Eq.lO, STDDM reduces to the second RPA (SRPA) 6]. The strength 
function S{E) defined as 

SiE)^ J2 \('f,mo)\'SiE~u;,) (4) 

for an excitation operator Q is given in terms of the solution of Eq.|(21). The detailed expressions of S{E) for one-body 
and two-body excitation operators are given in Ref.^5]. 

The Lipkin model 7] describes an N-fermions system with two N-fold degenerate levels with energies e/2 and — e/2, 
respectively. The upper and lower levels are labeled by quantum number p and — p, respectively, with p = 1^2, N. 
The model is described by the following hamiltonian 

H^eJ, + ^{Jl + J^), (5) 

where the operators are given as 

1 ^ 

J;, = - ^(a+flp - a_p+a_p) (6) 
p=i 

N 

j+ = ji^Y.4''-p- (7) 

p=i 

Using an iV = 4 system as an example, we solve Eqs.Q and Q for the ground states and Eq.® for the excited 
states. The hamiltonian matrix in STDDM (Eq.lO) is not hermitian in general 0,0 and, as a result, there is a 
possibility that some eigenvalues become complex. However, the hamiltonian in the Lipkin model is so simple that 
all eigenvalues are real and that the strength functions are positive definite. 

First we present the results for the ground states. We fix e at 2 and change V. To obtain converged results in 
the iterative gradient method, we need to start with a small value of interaction strength V and gradually increase 
it: We start with V/50 and increase it by V/50 for each 200 iterations. We found that the correlated ground states 
obtained coincide with those calculated using a time-dependent approach which provides us with another method 
for obtaining a stationary solution of the TDDM equations. The total energies obtained in TDDM (dashed line) for 
various interaction strength \V\ are shown in Fig. ^in comparison with the exact (solid line) and HF (dotted line) 
ones. The Lipkin model doesn't have stable HF ground states for \V\/e > 1/(N — 1). For \V\/e > 1/3 in the case of 

= 4, the 'deformed' HF solution is shown The results in TDDM and STDDM are always calculated using 

the original (not deformed) single-particle basis. The results in TDDM are very close to the exact ones in the wide 
range of the interaction strength: The deviation of the total energy from the exact one is only 1.1% at \V\/€ = 0.3. 

Now we show the results for the first excited state (one-phonon state). We use V = —0.6 which is slightly smaller 
than the critical value V — —2/3. Fig. [21 shows the strength function for the excitation operator Qi = J+ + J_ 
calculated in STDDM (dashed line). The strength functions for the exact solution (solid line) and the RPA solution 
(dotted line) are also drawn for comparison. The dot-dashed line depicts the result in a modified STDDM (mSTDDM) 
which will be discussed in connection with the two-phonon state. To make differences among various calculations easy 
to see, we smoothed the strength functions with an artificial width P/e = 0.025. RPA overestimates the collectivity 
of the first excited state: The excitation energy is much lower than the exact solution and the transition strength 
becomes quite large. On the other hand both the excitation energy and the transition strength calculated in STDDM 
are very close to the exact values. There appears a state at E/e = 1.3 in the STDDM result. This is due to the 
coupling to 3p-lh and lp-3h configurations. The value of the energy-weighted sum mi of the transition strength is 
14.2 in STDDM while the exact value is 14.6. The value of mi in RPA is 15.2. The results in STDDM demonstrate the 
importance of ground-state correlations and the coupling to the two-body amplitudes which are completely missing 
in RPA. The eff'ects of ground-state correlations may be classified as the self-energy contributions, the modification 
of p-h interactions and the vertex corrections For example, the fact that the excitation energy in STDDM is 

much larger than that in RPA is due to the increase in the self-energy of Ip-lh configurations. However, it is not 
easy in our approach to clearly distinguish the various contributions of the ground-state correlations. The excitation 
energy of the one phonon state is shown in FigOl as a function of the interaction strength. It is clearly seen that 
STDDM avoids the breakdown of RPA in the region of the phase transition (|V^|/e > 1/3) by including the effect of 
ground-state correlations. 
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FIG. 1: Total energy Ejt as a function of \V\ for iV = 4. The solid, dashed, and dotted lines depict the exact solutions, the 
TDDM results, and the HF results, respectively. 
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FIG. 2: Strength distributions of the one-phonon state calculated in STDDM (dashed Hue) and RPA (dotted line). The exact 
solution is shown with the solid line. The strength functions are smoothed with an artificial width F/e = 0.025. 
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FIG. 3: Excitation energy of the first excited state as a function of the interaction strength. The dashed and dotted lines 
connect the results in TDDM and RPA, respectively, while the solid line indicates the exact solutions. 
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FIG. 4: Strength distribution of the two-phonon state of the Lipkin model calculated in STDDM (dashed line). The solid 
and dotted lines depict the exact solution and the result in SRPA, respectively. The result in mSTDDM is shown with the 
dot-dashed line. The strength functions are smoothed with an artificial width F/e = 0.025. 



The strength function in STDDM (dashed line) for the two-phonon state excited with Q2 — (J+ + J )^ is shown 
in Fig. 0] where the strength function of the exact solution (sohd line) and that in SRPA (dotted line) are also 
drawn for comparison. The interaction strength used isV = —0.6. In SRPA where only 2p-2h and 2h-2p amplitudes 
are considered, neither coupling to the one-body amplitudes nor interaction among 2p-2h configurations exist in the 
Lipkin model. Therefore, the result in SRPA is equivalent to the unperturbed one. The state in STDDM has much 
larger transition strength than SRPA. This is due to the fact that Xppipiipm , Xphp'h' and Xhh'h"h"' amplitudes in 
addition to Xpp'hh' and Xhh'pp' are taken into account in STDDM. However, the transition strength is overestimated 
and the excitation energy is lower than the exact solution. This is similar to the relation between the RPA result 
and the exact one for the one-phonon state. This means that the effects of ground-state correlations are not fully 
taken into account in STDDM for the two-phonon state. In fact, the self-energy terms for 2p-2h configurations are not 
included in STDDM because they are kinds of three-body effects. In Ref.Q we presented a prescription for taking into 
account the three-body effects in which missing terms are added to the hamiltonian matrix acting on the two-body 
amplitudes. (See Eq.(25) of Ref.Q for a detailed explanation of the prescription.) We call this version of STDDM 
a modified STDDM (mSTDDM). The result in mSTDDM is shown in Fig. Hwith the dot-dashed Unc. It becomes 
closer to the exact solution: The transition strength is reduced and the excitation energy is increased. Due to the 
coupling to the two phonon state, the first excited state is slightly modified as well, as shown in FigI21 

In summary, we applied the TDDM plus STDDM approach to the Lipkin model. It was found that this approach 
gives much better results for the ground-state energies and the first excited states (the one-phonon state) than HF and 
RPA. However, the second excited state (two-phonon state) calculated in STDDM was not so good as the one-phonon 
state. This is due to the lack of some ground-state correlations coming from 3-body effects. It was shown that the 
modified STDDM which includes the three-body effects can improve the result for the two-phonon state. 
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